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SIMULATION AND VISUALIZATION OF FACE SEAL MOTION STABILITY 

8Y MEANS OF COMPUTER 6ENERATE0 MOVIES 

by 1. Etsion and B. M. Auer 

National Aeronautics and Space Administration 
Lewis Research Center 
Cleveland, Ohio 44135 

A computer aided design method for mechanical face seals is described. 
Based on computer simulation, the actual motion of the flexibly mounted ele- 
ment of the seal can be visualized. This is achieved by solving the equa- 
tions of motion of this element, calculating the displacements in its vari- 
ous degrees of freedom vs. time, and displaying the transient behavior in 
the form of a motion picture. Incorporating such method in the design phase 
allows one to detect instabilities and to correct undesirable behavior of 
the seal. A theoretical background is presentea. Details of the motion 
display technique are described, and the usefulness of the method is demon- 
strated by an example of a noncontacting conical face seal. 

NOMENCLATURE 

C Q equilibrium center-line clearance 

D simulated center-line clearance 

F dimensionless axial force 

1 dimensionless moment of inertia 

K* support stiffness 

M dimensionless moment 

m dimensionless mass 

m* ring mass 

R simulated seal radius 

r^ ring radius of gyration 


r m seal mean radius 

r 0 seal outer radius 

t dimensionless time 

Z* axial displacement 

Z dimensionless displacement, Z*/C Q 

a* tilt angle 

o normalized tilt, ** r Q IC Q 

e* coning angle 

$ normalized coning, &* r 0 /C Q 

y* nutation 

y normalized nutation, y* r Q IC Q 

y 1 simulated nutation 

^ precession 

(d shaft angular velocity 

INTRODUCTION 

A mechanical face seal consists of two discs one of which is attached 
to the shaft and rotates with it while the other one is stationary, usually 
one of the two discs is flexibly supported to allow self alignment and 
tracking of the mating disc. In figure 1 a face seal with fixed rotor and 
flexibly mounted stator is shown. The stator can move axially and tilt 
about two orthogonal diameters while its circumferential rotation is pre- 
vented by anti-rotation locks, it can also move radially in two perpendicu- 
lar directions. Thus, the flexibly mounted stator has 5 degrees of free- 
dom. The motion of the flexibly mounted element is controlled by the forces 
acting upon it as well as its mass inertia, ano stiffness and damping prop- 
erties of the dynamic system. This motion can be quite complex and diffi- 
cult to visualize from simple, two dimensional plots of time variations in 
the various degrees of freedom. 


2 


A well designed seal is one in which the motion of the flexibly mounted 
element is stable. Stable motion means steady tracking of the fixed element 
by the flexibly mounted one. The optimum design is one that results in per- 
fect tracking where the two discs remain parallel with a constant separation 
at all times. As was mentioned above many factors affect the dynamic be- 
haviour of the seal and hence the motion of the flexibly mounted element. 

The sal can be either stable or unstable depending on the compatibility of 
its operating conditions and design parameters. Indeed, various sources of 
seal instabilities were observed and reported in the literature, for exam- 
ple, references [1-6]. 

Many seal failures could be avoided if the designer had an easy and 
fast means by which his design could be checked. It is the objective of 
this paper to describe a method, based on computer simulation, by which the 
motion of the flexibly mounted element of the seal can be visualized. A 
noncontacting conical face seal was selected for the purpose of demonstrat- 
ing the method. Nevertheness, this method is suitable for any other con- 
figuration as long as the fluid film pressure distribution and the dynamic 
properties of the seal system are known. Basically the computer simulation 
is based on a solution of the equations of motion of the flexibly mounted 
element. Accelerations, velocities and displacements in the various degrees 
of freedom are calculated and the resultant transient behaviour of the seal 
is displayed in the form of a motion picture. Observation of the seal 
motion allows one to detect instabilities or any other undesirable be- 
haviour. It further allows the designer to correct such undesirable situa- 
tions by changing his design parameters and observing the resulting effect 
on the seal's dynamic behavior. 
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THEORETICAL BACKGROUND 


The seal model is shown in figure 2. The seal seat is the fixed rotor 
rotating at a speed u about the z axis of an inertial reference system 
xyz. The seal ring is the flexibly mounted stator. It can move axially 
along the z axis and tilt about the x and y axes of the inertial 
reference system xyz. These motions are the three major degrees of freedom 
of the seal ring. Usually the ring has two more degrees of freedom consist- 
ing of radial displacements along the x and y axes. However, these are 
much more restricted than the previous three major displacements and will 
not be treated here. Thus, the ring displacements are 2* along the z 
axis, and a* and o* about the x and y axes, respectively. The re- 
sultant orientation of the seal ring in the inertial reference xyz can be 
described by a rotating coordinate system 123. In this system axis 3 is 
perpendicular to the plane of the stator, axis 1 remains in the xy plane, 
and axis 2 passes through the point of maximum clearance. Thus, the 
coordinate system 123 has a nutation y* and a precession $ as shown in 
figure 2. It should be noted that while the ring itself is prevented from 
any circumferential rotation by the seal's anti-rotation locks, the system 
123 is free to rotate with respect to the ring while the plane 12 remains 
always in the plane of the ring. 

An axial force F* and tilting moments M* and M* are acting upon 
the stator causing its motion in the three major degrees of freedom. The 
system of force and moments is contriouted by the pressure in the fluid film 
between stator and rotor, by the balance pressure on the back of the stator, 
and by various elements of the flexible support, for example, mechanical 
springs and secondary seal. Reference l 6] describes in detail the system of 
forces and moments for the model of a coned seal with perfectly aligned 
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rotor. The motion of this coned model will be shown In the following 
discussion but, as was mentioned in the introduction, any other seal model 
can be treated along the same line of approach. 

The general dimensionless equations of motion of the flexibly mounted 
stator are 

F z « mZ (1) 

m x . is x 12) 

- I.y 13) 

Knowing the force F and moments M and M for a given position 

4 A J 

of the stator relative to the rotor, enables one to calculate the accelera- 
tions Z, o , and cl from equations (1) to (3). A time integration 
routine can then be used to find the velocities Z, a , a , and the new 

a jr 

displacements Z, <* , and a v . From the new relative position between 

a Jr 

stator and rotor a new system of force F z and moments M x and M y 
is found and the whole process repeats itself. 

The pressure contributed by the fluid film between stator and rotor is 
found by solving the Reynolds equation taking into account hydrostatic, 
hydrodynamic, and squeeze film effects [6]. Cavitation in the fluid film is 
also considered by summing up the hydrodynamic, squeeze, and hydrostatic 
pressures at each point in the sealing area and checking the total against a 
selected value of a cavitation pressure. If the total pressure at a point 
is less than the cavitation value it is replaced by the selected cavitation 
pressure prior to integration. 

Since the actual clearance in mechanical face seals is very small (of 
the order of few micrometers) the angles a*, a* and y* are also very 
small and, hence, can be treated as vector-. From figure 3 we have 

o* » Y* cos 1 J 1 (4) 
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y* sin ¥ 



( 5 ) 


Hence, the dimensionless nutation y and the precession angle ^ can Oe 
found from equations (4) and (5) in the form 



( 6 ) 

( 7 ) 


These two time dependent angles along with the axial oisplacement Z allow a 
three dimensional representation of the stator motion. Using computerized 
graphics this motion can be either watched on-line, or filmed frame by frame 
at each time step to provide a computer generated movie. 

SEAL MOTION OISPLAY 

The rotating seal seat and the flexibly mounted seal ring are re- 
presented in the movies by two circles, one above the other. The lower cir- 
cle represents the rotating seat and the upper circle represents the ring. 
The inertial coordinate system xyz with its origin at the center of the 
lower circle, as shown in figure 4, is usea as a reference to describe the 
relative motion between the seal seat and the seal ring. The rotation of 
the seat is represented by a moving symbol which is drawn on the lower cir- 
cle. This symbol represents the position of a fixed point on the seat as 
the seat rotates at a constant shaft speed. 

Both circles are drawn with an aroitrary radius k representing the 
outer radius of the seal. The origins of the two circles, designated by two 
different symbols, are located on the z axis. The design clearance is 
represented by an arbitrary distance 0 between the fixed origin of the 
lower circle and another fixed symbol on the z axis. 
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The equations for the bottom circle In the xyz coordinate system are 
the equations of a circle, namely, 

x . R cos d 
y m R sin i 

2 m 0 

where 0<: d < 2w and Is measured from the x axis. The equations for the 
upper circle, representing the seal ring, are more complex and are obtained 
by transformations through the precession angle and tne nutation angle 
y*. The procedure is as follows: 

Consider first a rotation of the xyz system of figure 2 by an angle 
i|; about the z axis. The new system designated xyz has its x and 2 
axes coincide with axes 1 and z, respectively, of figure 2. Any point in 
the xyz coordinate system can oe transformed into the xyz system by 


X * X COS 

♦ - y sin ^ 

(«) 

y * x sin 

'l' + y cos d 

(y) 

2 

* z 

(10) 


The coordinate system xyi is now rotated by an angle y* about the x 
axis. The new system is designated x*y*z* and its axes coincide with the 
coordinate system 123 of figure 2. The transformation from the coordinate 
system x*y*z* into the system xyz is 

7 = x* UD 

y * y* COS y* - Z* Sin y* ( 12 ) 

z « y* sin y* + z* cos y* (13) 
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Tne equations jf a circle in tne coordinate x*ste* are 

x* • R cos • 

>• - R sin # 

2 * • 0 

However, since tne seal ring Is prevented from rotation, ano tne coordinate 
system x*y*z* nas been rotated by an angle , with respect to the met 
tlal system xyz, any point on the upper circle has to be rotated ba*x l»> an 
angle (-\p) aoout tne z* axis. Thus, each point tK,#) on the upper t ir» le 
will nave, after the rotation, the coordinates 


R cos 4 cos <P ♦ R sin 4 sin 4 

(14) 

-R cos 4 sin <P ♦ R sin 4 cos <p 

(16) 

2* - 0 

(It) 


Using equations (11) to (16) in equations (Uj to (10), recalling that the 
origin of the upper circle in figure 4 is a distance 0(1 ♦ /) above Hie 

origin of the lower circle, the equations of the upper circle in the ay/ 

coordinate system of figure 4 are 
x - (R cos 4 cos <P ♦ R sin 4 sin <P)cos <P 

♦ (R cos 4 sin <p - H sin 4 cos <p)cos »' sin 4 
y • (R cos 4 cos <P ♦ R sin 4 sin <P)s1ri 4 

- (R cos 4 sin 4 - K sin 4 cos igcos r' cos 4 

z ■ 0(1 * l) - (R cos 4 sin 1» - R sin 4 cos *isln »' 

where y' It a scaled angle of nutation related to y* 0/ 

k sir. ,• r o lln '* 

0 ' l 0 

Since tne angle y* is ter/ small «e nare sin y* - ,* and 0 / using to* 
oimersionless angle of nutation, y, of equation (6>, eq^atlor '/'t, w 

T ‘ • Sin -1 (y jf) </, t 


nn 

( l hi 

n;, 

> /<), 


Because the physical angle y* and the design clearance C 0 are 
very small, 0 and R are chosen so that differences In the axial separa- 
tion and In the angle of nutation appear In the movie to oe relatively 
large. Thus, making it easier to visualize the seal motion. Care must be 
taken to choose R and U so that yD/R Is not greater than one. Typical 
values used were R « 2, 0 ■ 1. 

The x, y, and z coordinates defining the upper and lower circles of 
figure 4 have to oe transformed to a n coordinate system so that the 
circles can be projected onto a two dimensional frame as shown in figure b. 
The computer program uses the three spherical coordinate parameters R^, 

Op, and 0p of figure 5 to define the three dimensional perspective of 
the two circles. The circles are then viewed from a point in space defined 
by these three parameters where kp is the distance from the origin of 
the xyz coordinate system to the viewer's eye. In the present movies 
Rp was chosen to be 15, or 7.5 times the radius of the circle. Tne 
angles Op and dp were chosen to be both 45°. 

The displacements o , and a corresponding to times 

i y i 

t^, where t^ « t^ + at, i * 1, 2, ... n, were obtained from a solu- 
tion of equations (1) to (3J. The values of times and corresponding dis- 
placements were stored on a tape which was later used for the plotting pro- 
gram. Various values of the time step at were examined, and it was found 
that fifty time steps per one shaft revolution clearly reveal the seal 
motion. A varying number of points was tried to define the circumference of 
each circle. It was found that 40 circumferential points connected by 
straight lines are enough to obtain smooth circles. 
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The movie plotting program was run on an ! 3M 360 computer, and it took 
10 minutes of computer time to generate 1000 frames of film which correspond 
to 20 shaft revolutions. Each frame appeared on the display screen for 0.09 
seconds. 

Three examples consisting of portions of three computer generated 
movies are shown in figures 6 to 8. The frame sequence shown in these exam- 
pies is 5 frames per one shaft revolution. The number of shaft revolutions 
counted from the instance of an initial disturbance given to the upper cir- 
cle is shown on the lower left corner of each frame. The three examples 
correspond to the three basic modes of operation that were found in refer- 
ence [6] namely, unstable, critically stable and stable modes. In all 
three cases shown in figures 6 to 8 the seal radius ratio is r^/r Q - 0.9 
and the coning is **r Q / C Q « 10. All three cases have the same pressure 
parameter li‘ m /r sp ) 2 (p Q - Pj)r 2 /K*C 0 « 90.25, but they differ in 
their speed parameter values. The speed parameter values ( r g / r S p) 2 m*« 2 /K*, 
are 32, 21.95, and 8 for the cases shown in figures 6 to 8, respectively. 
Figure 9, which is taken from reference j.6 j shows stability maps for seals 
of radius ratio r^/r Q « 0.9. From this figure it is clear that the 
seal of figure 6 is unstable, the seal of figure 7 is critically stable, and 
the seal of figure 8 is stable. Indeed, in figure 6 a small initial dis- 
turbance is developed into a catastrophic failure, and rubbing contact oc- 
curs after about 20 shaft revolutions. The development of stator motion is 
shown in figure 6 for the 1st, 17th, 18th, and 19th revolutions. At the 
beginning the nutation angle increases very slowly but towards the 18th re- 
volution this increase becomes quite significant and is followed by a sub- 
stantial increase in the center- line clearance. 
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In the example of the critically stable seal shown In figure 7, the 
initial disturbance results in some angular and axial vibrations which are 
damped very rapialy. After about two shaft revolutions the stator assumes a 
constant nutation angle and from there on it wobbles at a constant frequency 
of half the shaft speed. Although the two seal faces do not contact, and 
hence no rubbing-caused failure occurs, the seal may still fail due to an 
excessive leakage caused by the large relative tilt between its mating 
faces, and the increase in the center-line clearance which are visible in 
figure 7. 

Figure a shows an example of a stable seal. Here a fairly large ini- 
tial disturbance causes axial and angular vibrations which are damped out 
very rapidly. After about two shaft revolutions the seal faces became 
parallel and the designed center-line clearance is maintained. 

CUNCLUDlhb REMARKS 

The simulation and visualization of noncontacting face seal motion 
stability by means of computer generated movies is described. Tne technique 
is based on solving the equations of motion of the flexibly mounted seal 
ring ano displaying the relative position of the seal ring and seal seat at 
subsequent time steps. Using computer graphics a perspective view of the 
seal is obtained on a display screen thus enabling observation of the seal 
motion. 

This simulation and visualization technique allows the designer to de- 
tect instabilities or any other undesirable behaviour of the seal. Thus, 
the designer can correct undesirable situations by changing the design 
parameters and observing their effect on the seal dynamics. It is believed 
that using such techniques can save time and money spent on testing, and can 
result in more reliable seal designs. 
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Figure 1. - Schematic of a radial face seal. 
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Figure 2. - Seal model and coordinates systems. 
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Figure 7. - Simulated seal motion at indicated no. of revolutions, critically stable case. 


[/ T~ ) ; 

A \ ■ J \ ■ j 


0.0 


0.2 

0.4 


0.6 

0.8 


v\ 

v . 

• \ 


\ .«* 

\ 

1.0 


12 

1.4 


L 6 

1.8 

• 

J 

v 

V. 

1 

•* 


2.0 


Z2 

Z4 


2.6 

2.8 

: 

• 

i * 

j [ 

V 

V 

V 

3.0 


3.2 

3.4 


3.6 

3.8 




oi* u. :i 


i»m;e is 

. WLITX 


Figure 8. - Simulated seal motion at indicated no. of revolutions, stable case. 




1 Report No 

NASA TM-81581 


2 Government Accession No 


3 Recipient % Catalog No 


4 Title and Subtitle 

SIMULATION AND VISUALIZATION OF FACE SEAL MOTION 
STABILITY BY MEANS OF COMPUTER GENERATED MOVIES 


7 Authotlsi 

I. Etsion and D. M. Auer 


9 Performing Organization Name and Address 

National Aeronautics and Space Administration 
Lewis Research Center 
Cleveland, Ohio 44135 


12 Sponsoring Agency Name and Address 

National Aeronautics and Space Administration 
Washington, D.C. 20546 


6 Refiort Date 

6 Performing Organi/jPon Code 


8 Performing Organization Report No 

E -554 


10 Work Unit No 


1 1 Contract or Grant No 


13 Type of Report and Period Covered ! 

Technical Memorandum 


14 Sponsoring Agency Code 


15 Supplementary Notes 

Prepared for the Fluid Sealing Conference sponsored by the British Hydrodynamic Research 


Association, Leeuwenhorst, Holland, April 1-3, 1981. 

16 Abstract 

A computer aided design method for mechanical face seals is described. Based on computer 
simulation, the actual motion of the flexibly mounted element of the seal can be visualized. 

This is achieved by solving the equations of motion of this element, calculating the displace- 
ments in its various degrees of freedom vs. time, and displaying the transient behavior in the 
form of a motion picture. Incorporating such method in the design phase allows one to detect 
instabilities and to correct undesirable behavior of the seal. A theoretical background is pre- 
sented. Details of the motion display technique are described, and the usefulness of the method 
is demonstrated by an example of a noncontacting conical face seal. 


17 Key Word' (Suggested by AuthorlsM 

Seals 

Mechanical face seals 
Seal dynamics 

18 Distribution Statement 

Unclassified - unlimited 
STAR Category 37 

19 Security Classif (of this report! 

Unclassified 

20 Security Clasuf (of this page' 
Unclassified 

21 No of Pages 

22 Price* 


* For sale by the National Technical Information Service Springfield V : g < n i a 22161 














